home *** CD-ROM | disk | FTP | other *** search
/ Amiga Format CD 46 / Amiga Format CD46 (1999-10-20)(Future Publishing)(GB)[!][issue 1999-12].iso / -in_the_mag- / reader_requests / scilab / tests / pow.dia.ref < prev    next >
Text File  |  1999-09-16  |  3KB  |  147 lines
  1.  
  2. //test with scalar vector and matrices
  3.  
  4. x1=1:3;
  5.  
  6. x0=0:2;
  7.  
  8. if or(x1**2<>[1 4 9]) then bugmes();quit;end
  9.  
  10. if or(x0**2<>[0 1 4]) then bugmes();quit;end
  11.  
  12. if or(x1**0<>[1 1 1]) then bugmes();quit;end
  13.  
  14.  
  15. if norm(x1**(2*(1+%eps))-x1**2)>100*%eps then bugmes();quit;end
  16.  
  17. if norm(x0**(2*(1+%eps))-x0**2)>100*%eps then bugmes();quit;end
  18.  
  19.  
  20. if norm(x1**(-2)-[1 0.25 1/9])>100*%eps then bugmes();quit;end
  21.  
  22. if norm(x1**(-2*(1+%eps))-[1 0.25 1/9])>100*%eps then bugmes();quit;end
  23.  
  24.  
  25. p=2+%eps*%i;
  26.  
  27. if norm(x1**p-x1**2)>100*%eps then bugmes();quit;end
  28.  
  29. if norm(x0**p-x0**2)>100*%eps then bugmes();quit;end
  30.  
  31. if norm(x1**(-p)-x1**(-2))>100*%eps then bugmes();quit;end
  32.  
  33.  
  34. y=%eps*%eps*ones(1,3);x1=x1+y;x0=x0+y;
  35.  
  36. if norm(x1**2-[1 4 9])>100*%eps then bugmes();quit;end
  37.  
  38. if norm(x0**2-[0 1 4])>100*%eps then bugmes();quit;end
  39.  
  40.  
  41. x1**2.000000001;x0**2.000000001;
  42.  
  43. if norm(x1**2-[1 4 9])>100*%eps then bugmes();quit;end
  44.  
  45. if norm(x0**2-[0 1 4])>100*%eps then bugmes();quit;end
  46.  
  47.  
  48. if norm(x1**(-2)-[1 0.25 1/9])>100*%eps then bugmes();quit;end
  49.  
  50. if norm(x1**(-2*(1+%eps))-[1 0.25 1/9])>100*%eps then bugmes();quit;end
  51.  
  52. if norm(x1**p-x1**2)>100*%eps then bugmes();quit;end
  53.  
  54. if norm(x0**p-x0**2)>100*%eps then bugmes();quit;end
  55.  
  56. if norm(x1**(-p)-x1**(-2))>100*%eps then bugmes();quit;end
  57.  
  58.  
  59. x1=[1 2;3 4];
  60.  
  61. if or(x1**1<>x1) then bugmes();quit;end
  62.  
  63. if or(x1**(-1)<>inv(x1)) then bugmes();quit;end
  64.  
  65. if or(x1**2<>x1*x1) then bugmes();quit;end
  66.  
  67. if or(x1**(-2)<>inv(x1)**2) then bugmes();quit;end
  68.  
  69.  
  70. x1(1,1)=x1(1,1)+%i;
  71.  
  72. if or(x1**2<>x1*x1) then bugmes();quit;end
  73.  
  74. if or(x1**(-2)<>inv(x1)**2) then bugmes();quit;end
  75.  
  76. if or(x1**1<>x1) then bugmes();quit;end
  77.  
  78. if or(x1**(-1)<>inv(x1)) then bugmes();quit;end
  79.  
  80.  
  81. if or(rand(4,4)**0<>eye(4,4)) then bugmes();quit;end
  82.  
  83.  
  84.  
  85. x1=[1 2;2 3];x2=x1**(1/2);
  86.  
  87. if norm(x2**(2)-x1)>100*%eps then bugmes();quit;end
  88.  
  89. x2=x1**(-1/2);
  90.  
  91. if norm(x2**(-2)-x1)>100*%eps then bugmes();quit;end
  92.  
  93.  
  94. //x1=[1 2+%i;2-%i 3];
  95.  
  96. //x2=x1**(1/2);
  97.  
  98. //if norm(x2**(2)-x1)>100*%eps then bugmes();quit;end
  99.  
  100. //x2=x1**(-1/2);
  101.  
  102. //if norm(x2**(-2)-x1)>100*%eps then bugmes();quit;end
  103.  
  104.  
  105. //test with polynomial vector and matrices
  106.  
  107. s=poly(0,'s');
  108.  
  109. if or(coeff(s**3+1)<>[1 0 0 1]) then bugmes();quit;end
  110.  
  111.  
  112.  
  113. x1=[1 s+1 s**3+1];
  114.  
  115. if or(x1**2<>[1 1+2*s+s**2  1+2*s**3+s**6]) then bugmes();quit;end
  116.  
  117. if or(coeff(x1**0)<>[1 1 1]) then bugmes();quit;end
  118.  
  119. if or(x1**3<>[1,1+3*s+3*s^2+s^3,1+3*s^3+3*s^6+s^9]) then bugmes();quit;end
  120.  
  121. if or((x1**(-1)-[1 1/(1+s)  1/(1+s**3)])<>0) then bugmes();quit;end
  122.  
  123.  
  124. x1=[s+1 2*s;3+4*s^2 4];
  125.  
  126. if or(x1**1<>x1) then bugmes();quit;end
  127.  
  128. if or(x1**(-1)<>inv(x1)) then bugmes();quit;end
  129.  
  130. if or(x1**2<>x1*x1) then bugmes();quit;end
  131.  
  132. if or(x1**(-2)<>inv(x1)**2) then bugmes();quit;end
  133.  
  134.  
  135. x1(1,1)=x1(1,1)+%i;
  136.  
  137. if or(x1**2<>x1*x1) then bugmes();quit;end
  138.  
  139. //if or(x1**(-2)<>inv(x1)**2) then bugmes();quit;end  //simp complexe non implemented
  140.  
  141. if or(x1**1<>x1) then bugmes();quit;end
  142.  
  143. //if or(x1**(-1)<>inv(x1)) then bugmes();quit;end //simp complexe non implemented
  144.  
  145.  
  146.  
  147.